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The (Ca a; Lai_ :E )(Bai.75_ 2 :Lao.25+i)Cu30 y system is ideal for testing theories of high temperature 
superconductivity, since nearly the full range of doping is controlled by y, and T™ ax is continuously 
controlled by x, with minimal structural changes. We investigate this system with both transverse 
and longitudinal field /j, SR. This allows us to re-examine the Uemura relation, the nature of the 
spontaneous magnetic fields below T c , and the relation between their appearance temperature T g and 
j.max Q ur ma j or findings are: (I) the Uemura relation is respected by CLBLCO more adequately 
than by other cuprates, (II) T g and T c are controlled by the same energy scale, (III) the phase 
separation between hole poor and hole rich regions is a microscopic one, and (IV) spontaneous 
magnetic fields appear gradually with no moment size evolution. 



The Uemura relation is a milestone in the research of 
high temperature superconductors [HTSC] and a key in=. 
gredient in most modern theories of HTSC. This relationtl 
states that in underdoped HTSC the superconductivity 
transition temperature T c is proportional to the muon 
spin rotation [/iSR] line width a, and that the proportion- 
ality constant is universal for all HTSC materials. Since 
a is proportional to where X a b is the in-plane pene- 
tration depth [see Sec. |j, there is a one to one relation 
between T c and A~ & 2 for all underdoped HTSC. However, 
this relation seems to break in the optimally and over 
doped regions. To date, the origin of this breakdown is 
not clear. 

Another key ingredient of modern theories of HTSC 
is that at low temperatures (T), cuprates phase-jSeparate 
into regions that are hole "poor" and hole "rich"0. While 
hole rich regions become superconducting below T c , the 
behavior of hole poor regions at these temperatures is not 
quite clear. Some data support the existence of magnetic 
moments in these regionso, but there are still many open 
questions regarding these moments and the spontaneous 
magnetic fields associated with them. For example: Is 
there a true phase transition at T g l What is the field pro- 
file and how is it different from, or similar to, a canonical 
spin glass? Is the field confined solely to the hole poor 
regions or does it penetrate the hole rich regions? Also, 
the interplay between magnetism and superconductivity 
is not clear. Is strong magnetic background beneficial or 
detrimental to superconductivity? 

We address these questions by investigating the 
(Ca x Lai_ x )(Bai.75_ x Lao.25+i)Cu 3 O a (CLBLCO) family 
of superconductors. These superconductors belong to 
the 1:2:3 family and has several properties that make 
it ideal for our purpose. It is tetragonal throughout its 
range of existence < x < 0.5, so thexe is no order- 
ing of CuO chains. Simple valency sumsa, more sophis- 
ticated bond- valance, calculations^, and thermoelectric 
power measurements^] show that the hole concentration 
is x independent. As shown in Fig. [I], by changing y, for 



a constant value of x, the full superconductivity curve, 
from the underdoped to the overdoped, can be obtained. 
Finally, for different Ca contents, parallel curves of T c 
vs y are generated. Therefore, with CLBLCO one can 
move continuously, and with minimal structural changes, 
from a superconductor resembling YBCO to one similar 
to LSCO. We study the superconducting and magnetic 
properties of CLBLCO by performing, zero, longitudinal, 
and transverse field muon spin relaxation experiments. 
The report on the zero and longitudinaLfield measure- 
ments is an extension of our recent lettenJ. 
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FIG. 1. Phase dia- 

gram of (Ca :c Lai- : t)(Bai.75-xLao.2B+a:)Cu30 J/ . The dashed 
lines indicate samples with equal T g . 



I. EXPERIMENTAL ASPECTS 

The preparation of the samples is described elsewhere!. 
Oxygen content was determined using iodometric titra- 
tion. All the samples were characterized using X-ray 
diffraction and were found to be single phase. T c pre- 
sented in the phase diagram is obtained from resistivity 
measurements. 
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FIG. 2. jj.SK experimental setup representing the (a) lon- 
gitudinal and transverse field configurations in the (b) ISIS 
and (c) PSI facilities. 

Our /xSR experiments were done at two facilities. 
When a good determination of the base line was needed 
we used the ISIS pulsed muon facility at Rutherford 
Appleton Laboratory, UK. When high timing resolu- 
tion was required we worked at the Paul Scherrer In- 
stitute, Switzerland (PSI). In /iSR experiments one in- 
jects polarized muons into the sample and applies a 
magnetic field. Decay positrons are emitted preferen- 
tially in the direction parallel to the muon spin and 
are detected by positron counters. We used two types 
of setups, the transverse field [TF] and the longitudi- 
nal field [LF], which also includes zero field, as shown 
in Fig. ^. The LF configuration works as in panel (a) 
in both facilities. The TF configurations work in ISIS 
as in panel (b) and at PSI as in panel (c). From the 
counters an histogram is generated of positrons counts 
as a function of time, where the time is measured from 
the moment the muons enter the sample. The asym- 
metries Apsit), AuE>(t) and A^^t) are than generated 
by A FB {t) = [F(t] - B{t)] I [F(t) + B(t)] and similarly 
with U, D instead of F, B, where F, B, U, D are the 
forward, backward, up, and down counters. ALn(t) is 
generated only from the left counter in Fig. ||(c) . Most of 
the data were taken with a 4 He cryostat. However, in or- 
der to study the internal field profile we had to avoid dy- 
namical fluctuations by freezing the moments completely. 



For this purpose we used the 3 He cryostat at ISIS with 
a base temperature of 350 mK. All fiSH measurements 
were done on sintered pellets. 

II. TF-/1SR 

These experiments are done by field cooling (FC) the 
sample to 1.8 K at an external field of 3 kOe in PSI and 
400 Oe in ISIS. As explained above we apply the field per- 
pendicular to their spin direction, and every muon then 
precesses according to the local field in its environment. 
When field cooling the sample, a vortex lattice is formed, 
and the field from these vortices decay on a length scale 
of A. This leads to a inhomogeneous field distribution 
in the sample. Since the magnetic length scale is much 
larger than the atomic one, the muons probe the mag- 
netic field distribution randomly, which, in turn, leads to 
a damping of the muons average spin polarization. This 
situation is demonstrated in Fig. |3J where we present an 
image of the field profile, and the corresponding the real 
and imaginary part of the muon asymmetry. At temper- 
atures above T c the field is homogeneous and all muons 
experience the same field, and therefore no relaxation is 
observed. Well below T c there are strong field variations 
and therefore different muon precess with different fre- 
quencies, and the average polarization quickly decays to 
zero. In intermediate temperatures the field variation are 
not severe and the relaxation is moderate. 
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FIG. 3. Demonstrating the relation between the field distri- 
bution and the real and imaginary asymmetries in a TF-/iSR 
experiment. This data was taken at ISIS and is presented in 
a rotating reference frame. 

It was shown that in powder samples gf HTSC the 
muon asymmetry A(t) is well described byO, 

A(t) = cxp(-cr 2 i 2 /2) cos(wi + tp) (1) 

where lu — ^y^H is the precession frequency of the muon, 
a is the relaxation rate, and tp is a phase which depends 
on the counters used to generate the asymmetry Our 
analysis for both ISIS and PSI data is done in a ref- 
erence frame rotating at u rr f and the real and imagi- 
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nary components of the signal are fitted simultaneously. 
Therefore, the frequency in Fig. || is 7 M i? — u> rr f where 
Lj rr f is chosen arbitrarily for presentation purpose. The 
solid line in this figure is the fit result. The fact that 
the whole asymmetry relaxes indicates that CLBLCO is 
a bulk superconductor. 

The fit results for a are shown in Fig |[ As can be 
seen, the dependence of T c on a is linear in the under- 
doped region and universal for all CLBLCO families, as 
expected from the Uemura relations. However, there is 
a new aspect in this plot. There is no "boomerang" ef- 
fect, namely, overdoped and underdoped samples with 
equal T c have the same tr, with only slight deviations for 
the x = 0.1 sample as demonstrated by the arrows in 
Fig. ^. Therefore, in CLBLCO there is one to one corre- 
spondence between T c and cr, and therefore A~ b 2 , over the 
whole doping range. This is our first important finding. 
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FIG. 4. A Uemura plot showing T c vs. the muon relaxation 
rate a in Eq. |l| for the CLBLCO family of superconductors. 

III. ZF-^SR 

Typical muon asymmetry depolarization curves are 
shown in Fig. |5| (a) for different temperatures in the 
x = 0.1 and y = 7.012 (T c = 33. IK) sample. The change 
of the polarization shape with temperature indicates a 
freezing process, and the data can be divided into three 
temperature regions. In region (I), given by T > 8 K, the 
muon relaxes according to the well known Kubo-Toyabe 
(KT) function given by 

KT{t) = 1 + 1(1- A 2 t 2 ) exp(-iA 2 t 2 ), (2) 

(see Eq. ^) typical ofj-the case where only frozen nuclear 
moments are presentfl. In region (II), bounded by 8 K 
>T>3K, part of the polarization relaxes fast and the 
rest relaxes as in the first region. As the temperature 
is lowered the fast portion increases at the expense of 
the slow one. Moreover, the relaxation rate in the fast 
portion seems independent of temperature. Finally, at 
long time the asymmetry relaxes to zero. In region (III), 
where 3 K > T, the asymmetry at long times no longer 



relaxes to zero, but instead recovers to a finite value. 
This value is ~ 1/3 of the initial asymmetry A z (0). 
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FIG. 5. (a) ZF-/iSR spectra obtained in a x = 0.1, 
y = 7.012 sample with T C =33.1K. The solid lines are fit to 
the data using Eq. ^[ the dashed line is a fit using the sim- 
ulation as described in the text, (b) pSR spectra obtained 
in longitudinal fields from the x = 0.4, y — 6.984 sample at 
350 mK. (c) Polarization curves generated by the simulation 
program as described in the text. 

To demonstrate that the internal field is static at base 
temperature, the muon polarization was measured with 
an external field applied parallel to the initial muon spin- 
polarization. This geometry allows one to distinguish 
between dynamic and static internal fields. In the dy- 
namic case the asymmetry is field independent^!! In con- 
trast, in the static case the total field experienced by the 
muon is a vector sum of H and the internal fields, which 

are of order (B 2 ) 1 ' 2 . For H > (B 2 ) 1 ' 2 the total field 
is nearly parallel to the polarization. Therefore, in the 
static case, as H increases, the depolarization decreases, 
and the asymmetry recovers to its initial value. Because 
we are dealing with a superconductor, this field sweep 
was done in field-cool conditions. Every time the field 
was changed the sample was warmed to above T c and 
cooled down in a new field. The results are shown in 
Fig. ||(b). At an external field of 250 G, the total asym- 
metry is nearly recovered. Considering the fact that the 
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internal field is smaller than the external one due to the 
Meissner effect, this recovery indicates that the internal 
field is static and of the order of tens of Gauss. Next 
we perform quantitative data analysis in two parts: high 
temperatures (region II), and base temperature. 



A. High T Data Analysis 

First we discuss region II. Here we focus on the deter- 
mination oiTg. For that purpose we fit a combination pfj 
a fast relaxing function and a KT function to the datatll 



A z (t) = A m exp (-Vxt) +A n KT(t), 



(3) 



where A m denotes the amplitude of the magnetic part, A 
is the relaxation rate of the magnetic part, and A n is the 
amplitude of the nuclear part. The relaxation rate of the 
KT part was determined at high temperatures and is as- 
sumed to be temperature independent. The sum A 7n -\- A n 
is constrained to be equal to the total initial asymmetry 
at high temperatures. The relaxation rate A is common 
to all temperatures. The solid lines in Fig. ||(a) are the 
fits to the data using Eq. ||. 
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FIG. 6. Magnetic amplitude as function of temperature f 
or different samples. The solid lines are guides to the eye. 

The success of this fit indicates the simultaneous pres- 
ence of two phases in the sample; part of the muons probe 
the magnetic phase while others probe only nuclear mo- 
ments. As the temperature decreases A m , which is pre- 
sented in Fig. H for three samples, grows at the expense 
of A n . At low temperatures A m saturates to the full 
muon asymmetry. A similar temperature dependence of 
A m is found in all samples. The origin of the magnetic 
phase is electronic moments that slow down and freeze in 
a random orientation. The fact that A is temperature in- 

1/2 

dependent means that in the magnetic phase 7 M (-B 2 ) , 
where 7^ is the muon gyromagnetic ratio, is tempera- 
ture independent. In other words, as the temperature 
is lowered, more and more parts of the sample become 
magnetic, but the moments in these parts saturate upon 
freezing. 
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FIG. 7. T g vs. y. The horizontal solid lines are the equal 
T g lines appearing in Fig. |l|. 

Our criterion for T g is the temperature at which A m is 
half of the total muon polarization as demonstrated by 
the vertical lines in Fig. || for the three different samples. 
The phase diagram that is shown in Fig. ^ represents T g 
for various samples differing in Ca and O contents. This 
diagram is systematic and rather smooth suggesting good 
control of sample preparation. As expected, for constant 
x, higher doping gives lower T g . 

We have singled out three groups of samples with a 
common T g = 11, 8 and 5 K as shown in Fig. ^ by the 
horizontal solid lines. These samples are represented in 
the phase diagram in Fig. [j] by the dotted lines. The 
phase diagram, containing both T g and T c , is the second 
main finding of this work. It provides clear evidence of 
the important role of the magnetic interactions in high 
temperature superconductivity as discussed in Sec. [v|. 



B. Low T Data Analysis 

We now turn to discuss the muon depolarization at 
base temperature. In this case all the muons experi- 
ence only a static magnetic field, as proven above. This 
allows one to reconstruct the internal field distribution 
out of the polarization curve. The polarization of a 
muon spin experiencing a unique field B is given by 
P z (t) = cos 2 (6>) + sin 2 (6») cos(7|B|t), where 9 is the angle 
between the field and the initial spin direction. When 
there is an isotropic distribution of fields, a 3D powder 
averaging leads to 



PS) 



p(\B\)cos(j\B\t)B 2 dB 



(4) 



where p(|-B|) is the distribution of |B|. Therefore, the po- 
larization is given by the Fourier transform of p(\B\)B 2 
and has a 1/3 base line. When the distribution of B is 
centered around zero field, p(\B\)B 2 is a function with a 
peak at (B) and a width A, and both these numbers are 
of the same order of magnitude [e.g. Fig||(b)]. Therefore 
we expect the polarization to have a damped oscillation 
and to recover to 1/3, a phenomenon known as the dip 
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[e.g. the inset in Fig || (b)]. Gaussian, Lptjentzian and 
even exponential random field distribution^, and, more 
importantly, all known canonical spin glasses, produce 
polarization curves that have a dip before the 1/3 re- 
covery. This is demonstrated in Fig. O. For a Gaussian 
distribution of width A we obtain Eq.^ which is demon- 
strated in panel (a). The cases of a canonical spin glass 
Feo.05TiS2, and an extremely underdoped CLBLCO are 
presented in panels (b) and (c). Furthermore, a dipless 
polarization curve that saturates to 1/3 cannot be ex- 
plained using dynamical arguments. Therefore, the most 
outstanding feature of the muon polarization curve at 
base temperature is the fact that no dip is present, al- 
though there is a 1/3 tail. This behavior was found in 
all of qw samples with T c > 7r-K, and also in Ca doped 
YBCOEJ and Li doped YBColii 



[panel (a)]. Here (B) is smaller than A, and, indeed, the 
associated polarization in the inset is dipless. Thus in 
order to fit the base temperature polarization curve we 
should look for p(\B\)B 2 with most of its weight around 
zero field. This means that diverges like 1/B 2 at 

|-B| — 0, namely, there is abnormally high number of low 
field sites. 

It also means that the phase separation is not a macro- 
scopic one. If it were, all muons in the field free part 
would probe only nuclear moment and their polarization 
curve should have a dip or at least its beginning as in 
the high temperature data. The same would apply for 
the total polarization curve, in contrast to observation. 
Thus, the superconducting and magnetic regions are in- 
tercalated on a microscopic scale (~ 20A)E_3. This is the 
third main finding of this work. 
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FIG. 8. (a) The internal field distribution extracted from 
the simulations for the case of correlation length £ = 3 lattice 
constants, maximum moment size of 0.06/is and magnetic 
moment concentration p = 15%. Inset: The muon spin polar- 
ization for that distribution, (b) The same as above for the 
case of p = 35%. 

The lack of the dip in P z {t) can tell much about the 
internal field distribution. It means that (B) is much 
smaller than A. In that case the oscillations will be 
over-damped and the polarization dipless! In Fig. || we 
show, in addition to the (B) ~ A case described above 
[panel (b)], a field distribution that peaks around zero 



FIG. 9. Demonstrating the expected muon spin plariza- 
tion function for (a) a Gaussina field distribution, (b) in a 
canonical spin glass Feo.osTiSa, and (c) extremely underdoped 
CLBLCO. 

The special internal field distribution, and the nature 
of the gradual freezing of the spins, can be explained by 
the intrinsic inhomogeneity of hole concentration. The 
part of the sample that is hole poor, and for that reason 
is "more" antiferromagnetic, will freeze, while the part 
which is hole rich will not freeze at all. The variation in 
the freezing temperature of different parts of the sample 



5 



can be explained by the distribution of sizes and hole 
concentration in these antiferromagnetic islandsc3. The 
large number of low field sites is a result of the fact that 
the magnetic field generated in the magnetic regions will 
penetrate into the hole rich regions but not completely. 

IV. NUMERICAL SIMULATION 

To improve our understanding of the muon polariza- 
tion, we performed simulations of a toy model aimed at 
reproducing the results described above. 

A 2D 100 x 100 square lattice is filled with two kinds 
of moments, nuclear and electronic. All the nuclear mo- 
ments are of the same size, they are frozen and they point 
in random directions. Of the electronic moments only a 
small fraction p is assumed to be frozen; they represent 
magnetic regions with uncompensated antiferromagnetic 
interactions. Since these regions may vary in size, the 
moments representing them are random, up to a max- 
imum size. The frozen electronic moments induce spin 
polarization in the other electronic moments surrounding 
them. Following the work of otherstSl, we use decaying 
staggered spin susceptibility which we take to be expo- 
nential, namely, 

x(r) = (-l) n " +n »«p(-r/0 (5) 

where r = n x ax + n y ay represents the position of the 
neighbor Cu sites, a is the lattice vector, and £ is the 
characteristic length scale. Because of this decay, at low 
frozen spin concentration, large parts of the lattice are 
practically field free (expect for nuclear moments). How- 
ever, the important point is that no clear distinction be- 
tween magnetic and field free (superconducting) regions 
exists. This situation is demonstrated in Fig. [HJ. 



FIG. 10. Demonstrating the numerical simulations. Two 
spins (long arrows) are placed on the lattice. They polarize 
the near by spins. The muon interacts with the spin by dipo- 
lar interaction. Nuclear moment (which participate in the 
simulation) are not shown. 



The muon polarization time evolution in this kind of 
field distribution is numerically simulated. The interac- 
tion between the muon and all the other moments is taken 
to be dipolar, and £ is taken to be 3 lattice constantsBEZI. 
The dashed line in Fig. | is a fit to the T = 350 mK 
data, which yield p — 15% and maximum moment size 
~ 0.06/is ■ As can be seen, the line fits the data very 
well. However, as expected, the fit is sensitive to p£ 2 
only, namely the effective area of the magnetic islands, 
so longer £ would have given smaller p. The field distri- 
butions and the polarization curve shown in Fig. |8| were 
actually generated using the simulation. In (a) the spin 
density is 15% while in (b) the density is 35%. In panel 
(c) of Fig. H we show the spin polarization for different 
hole concentration, varying from 0% to 35% with the 
same £ = 3. The resemblance between the simulation 
results as a function of p and the muon polarization as 
a function of temperature in panel (a) leads us to our 
fourth conclusion that the freezing process is mostly a 
growth in the total area of the frozen AF islands. 

V. DISCUSSION 

We now discuss the phase diagram presented_ki Fig. [I]. 
This diagram is consistent with recent theoriesEj of hole 
pair boson motion in an antiferromagnetic background. 
Those theories conclude that 

T c oc Jn s (6) 

where n s is the superconducting carriep-density, and J 
is the antiferromagnetic coupling energycJ. Let us define 
Aphi = Pu — Po where phi is the number of mobile holes 
and po is the number of mobile holes at optimum. 
We assume that 

n s (Ap M ) = i (p + Ap h i) (7) 

and write 

= K(x)Ay (8) 

where Ay = y — 7.15 is chemical doping measured from 
optimum, and if is a scaling parameter which relates Ay 
to ApM. Since there is a linear dependence between T g 
and chemical doping (Fig. [?]) we predict that 

T g cx J(l - c g n s ). (9) 

^From Eq. | T c max cx Jn s (0), therefore both T c /T c max 
and T 9 /T™ ax should be functions only of Ap hl . This 
is demonstrated in Fig. |ll](a). We find K{x) by mak- 
ing all T c /r™ ax collapse pato one curve resembling the 
curve of La2- a; Sr a ;Cui04Lil, where the exact doping is 
known. Using these values of K(x) we also plot T g /T™ ax 
as a function of Ap^i in Fig. |ll|(b). Again all data 
sets collapse onto a single line described by T g /T™ ax = 
-3.1(2)Ap w - 0.21(2), or 
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T g = 0.3T" ax (l - c g n s ) (10) 

with c g = 10.3 when converting back to Eq. ^. This 
indicates that the same single energy scale J controls 
both the superconducting and magnetic transitions, and 
provides the exact n s and T™ ax dependence of T g . 



Finally, and most importantly, the phase diagram con- 
taining both T c and T g leads us to believe that these 
temperatures are determined by the same energy scale. 
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FIG. 11. T c /T™ ax and T g /T™ ax as a function of K(x)Ay 
where Ay = y- 7.15, and K(x) is chosen so that all T c /T™ ax 
data sets collapse to a single curve, which resembles the LSCO 
curve. 

VI. CONCLUSIONS 

We are now in a position to address the questions pre- 
sented in the introduction. The Uemura relations is ad- 
equately respected for our HTSC "families" . We believe 
this is a result of the fact that there are no structural 
changes between the different "families". The appear- 
ance of spontaneous magnetic held in CLBLCO is a grad- 
ual process. As the temperature is lowered microscopic 
regions of frozen moments appear in the samples, and 
their area increases but the moments do not. In the 
ground state the held profile is very different from that 
of a canonical spin glass or any other standard magnet. 
It could only be generated by microscopic intercalation 
of an abnormal number of zero held regions with mag- 
netic regions without a clear distinction between the two. 
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